Mechanical, Electrical, and Magnetic Properties of Ni Nanocontacts 
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The dynamic deformation upon stretching of Ni nanowires as those formed with mechanicaUy 
controllable break junctions or with a scanning tunneling microscope is studied both experimentally 
and theoretically. Molecular dynamics simulations of the breaking process are performed. In addi- 
tion, and in order to compare with experiments, we also compute the transport properties in the 
last stages before failure using the first-principles implementation of Landauer's formalism included 
in our transport package ALACANT. 



I. INTRODUCTION 

In a foreseeable future, the functionality of elec- 
tronic devices will rely on the conduction properties of 
molecules or nanoscopic regions comprised of a surpris- 
ingly small number of atoms. Over the past 10 years, var- 
ious experimental groups have developed different tech- 
niques to connect two large metallic electrodes by just 
an atom or a chain of atoms^'^''^. These systems receive 
names such as atomic-size contacts or nanocontacts. Al- 
though they are not expected to be of any practical tech- 
nological application in the near future, these systems are 
an excellent test bed to learn about electrical transport 
at the atomic scale. 

While a large amount of experimental and theoretical 
work has been reported for many metals, a deep theoret- 
ical understanding is still lacking in the case of magnetic 
nanocontacts, which exhibit a very rich and complex 
behavior^^. Modeling their mechanical, electrical, and 
magnetic properties with accuracy is a challenge from 
which we expect to learn important lessons on our way 
toward reliable theoretical descriptions of more sophis- 
ticated systems of relevance in present and future spin- 
based devices. We present here a comparison between 
theoretical results of the mechanical, magnetic, and con- 
duction properties of Ni nanocontacts, and experiments 
carried out in our laboratory. 



II. SCANNING TUNNELING MICROSCOPE 
HISTOGRAMS 



For the experiments, we used a high-stability scanning 
tunneling microscope (STM) at low temperatures (4.2K) 
and under cryogenic vacuum conditions. For both tip 
and sample, we used Ni wire 99.99+% pure. The wire was 
cleaned by sonication in an acetone bath and scratched to 
remove contamination attached to the surface. The con- 
ductance properties of the contacts formed in between 
tip and sample were measured in a typical two-probe 
configuration. A constant bias voltage (typically 10-100 
mV) was applied between tip and sample, and the cur- 
rent was measured using a homemade current amplifier in 
the range of tens of microammeters. We recorded traces 
of conductance as a function of the relative tip-sample 
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FIG. 1: Experimental conductance histogram for Ni nanocon- 
tacs recorded at a bias voltage of 100 mV and a temperature 
of 4.2K where two low-conductance peaks are clearly visible. 



distance as the two electrodes were brought together and 
separated. In every trace, we made a deep tip-sample in- 
dentation in order to prevent the repetition of the same 
atomistic configurations and to assure the clearness of the 
contacts. Afterward, the traces were collected to build a 
conductance histogram such as the one shown in Fig. 1. 

Histograms, similar to the aforementioned, for the first 
stages of conduction in Ni nanocontacts have been stud- 
ied before^i^i^. There a broad peak around 1.6 Go, where 
Go — 2e'^/h is the quantum of conductance, has been re- 
ported as the first peak after vacuum tunneling. This 
peak is attributed to the cases in which the contact con- 
sists of a single atom. On the other hand, for values of 
conductance below the one-atom peak, we can notice a 
large amount of data coming from tunneling. This effect 
is stronger in Ni than in other metals, such as Au, since 
for Ni, some of the traces show a smooth transition from 
tunneling to contact without a jumpi. 

Here, we have studied in detail the lowest conductance 
peak after tunneling and noticed that, indeed, it is not a 
single broad peak, but the superposition of two at around 
1.2 Go and 1.5 Go . The position of these peaks slightly 
changes for different contacts, and this may be the rea- 
son why they have not always been clearly resolved. We 
have performed separated conductance histograms from 
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FIG. 2: Inset: experimental breaking histograms recorded at 
4.2K at different bias voltages. The peak above 1 Go can be 
resolved as two different peaks, marked as lower conductance 
peak (LC) and higher conductance peak (HC). The figure 
shows the dependence of the position of the two peaks with 
the applied bias voltage averaged over different experimental 
realizations. 

the traces for the cases of either forming or breaking the 
contacts. There we find a different ratio in the height of 
the peaks, being in the case of breaking traces the peak 
at 1.2 Go , in general, higher than the one at 1.5 Gq and 
vice versa in the case of making the contacts. Finally, we 
notice a dependence of the position of the peaks with the 
bias voltage with a variation of even 0.3 Gq in a voltage 
range of 300 mV, as shown in Fig. 2. 

To the best of our knowledge, the described features 
have not been previously reported for any material. Nor- 
mally, the peaks in the conductance histograms do not 
change as a function of the bias voltage^, and are equal 
for the cases of breaking or forming the contacts. In or- 
der to understand our observations, one should first try to 
identify all the possible atomic configurations that could 
lead to a conductance in between 1 Go and 2 Go in Ni 
nanocontacts, and next look for the reasons that make 
the value of conductance to be so dependent on the bias 
voltage. In the following discussion, we will address the 
first question. 



III. MECHANICAL PROPERTIES 

The dynamic deformation of Ni nanowires upon 
stretching until failure has been studied using molec- 
ular dynamics with empirical potentials. This type of 
modeling has provided significant information about the 
atomic scale processes occurring during deformation of 
nanowires^i^. In this particular work, we focus on the 
last stages before failure of Ni nanocontacts. The inter- 
atomic potential for Ni developed by Mishin et al.— was 
used in these calculations. This potential has been fit- 
ted to reproduce the stacking fault energy of Ni. In all 



FIG. 3: Minimum cross section as a function of elongation 
along the [1 0] direction for a case with 2645 atoms. Insert 
shows configurations at different stages of the deformation. 



calculations, deformation is achieved by displacing the 
outer two layers of atoms on each side of the simulation 
box a fixed distance every 1000 simulation steps, simi- 
lar to what is done by other authors^ ^. Two different 
deformation velocities were used, 1 and 10 m/sii. 

The dependence of different parameters on the defor- 
mation and, in particular, on the last stages before fail- 
ure has been studied. On one hand, we have compared 
the deformation of different system sizes, between 77 and 
2645 atoms, with initial cross sections between 1.5 and 
3.5ao , where oq is the lattice parameter, ao =3.52 A , 
for the case where tension is applied along the [1 0] di- 
rection. The dependence with crystallographic direction 
has also been studied for systems with similar number of 
atoms (between 610 and 658) and for directions [10 0], 
[1 1 1], [1 10], and [1 1 2]. All these calculations were 
performed at a fixed temperature of 4.2K by rescaling 
the velocities of all atoms. Finally, the dependence with 
temperature was also studied for the particular case of 
deformation along the [1 0] direction and a cross sec- 
tion of 2ao. 

The minimum cross section perpendicular to the ap- 
plied tension is computed every 1000 steps following the 
method developed by Bratkovsky et al.^«. This method 
allows us to compare the results of deformation along 
different crystallographic directions. Fig. 3 shows the 
minimum cross section as a function of elongation ob- 
tained for one particular case with 2645 atoms and de- 
formation along the [1 0] direction. The insets show 
several configurations during the deformation: the ini- 
tial configuration, an intermediate configuration that is 
particularly stable, and the final configuration, which, in 
this case, consists of a single atom, a monomer, connect- 
ing the two sides of the nanowire. The mechanisms for 
deformation at this scale have been studied in detail by 
other authors^. Consistent with their work, we observe 
the sliding of planes during deformation (such as in the 
intermediate inset in Fig. 3) that results in a narrowing 
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FIG. 4: Histogram of minimum cross sections from 100 independent simulations for crystallographic directions, (a) [1 0] and 
[1 1 2]. (b) [1 1 1] and [1 1 0]. 
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FIG. 5: Histogram of minimum cross sections from 25 inde- 
pendent simulations and four different temperatures. [1 0] 
crystallographic direction. 



of the wire. This results in preferential configurations, 
reflected in the plateaus observed in Fig. 3. 

The dependence of the deformation on the stretching 
direction has been studied for [1 0], [1 1 1], [1 1 0], 
and [1 1 2] directions. Static calculations to obtain the 
energy to fracture along these different directions show 
that the [1 1 2] direction has the lowest energy per unit 
surface, followed by the [1 0], [1 1 0], and [1 1 1] direc- 
tions. However, this behavior could be different during 
dynamic deformation. For each direction, calculations 
were repeated 100 times in order to gain some statis- 
tics and obtain a histogram of cross sections. Fig. 4 
shows the histograms obtained for all crystallographic 
directions with clear peaks at particular cross sections. 
These preferential cross sections do not depend on the 
system size. Histograms obtained from smaller systems 
result in peaks located at exactly the same positions. It 
is interesting to point out that the first four peaks ap- 
pear in all cases. These correspond to the smallest cross 
sections, consisting of less than three atoms across. 

These first peaks are also very stable with temperature. 



Fig. 5 shows the histograms obtained for four different 
temperatures, 4.2K as before and 300K, 610K, and 770K. 
These histograms were obtained from 25 independent cal- 
culations and for the [1 0] direction. Notice that from 
4.2K to room temperature, there is a strong reduction in 
the first peak. Therefore, the contact breaks very rapidly 
at high temperatures. On the contrary, structures with 
cross sections of two or three atoms seem to be very sta- 
ble with temperature. For wider structures, there is not 
a clear peak at high temperature as in the case of 4.2K, 
which seems to point to very different types of structures 
possible when temperature increases. The dependence 
with temperature of these structures has been studied 
previously by other authors^. 

In what follows, we focus on the final stage before fail- 
ure of these nanocontacts. Two structures have been 
identified: a monomer, where a single atom acts as a 
bridge between the two contacts, and a dimer,where two 
atoms aligned forming a bridge between the two sides of 
the wire. These two configurations are shown in the in- 
sets of Fig. 6, (a) being the monomer and (b) the dimer. 
In order to identify the dimer, we have calculated the 
total number of neighbors for each atom and the number 
of neighbors to the left of the atom position and to the 
right along the z-direction. In this manner, it is easy to 
identify a dimer since it will consist of two neighboring 
atoms each one with only one neighbor on one side, one 
atom to the left, the second one to the right. From all the 
cases computed including all crystallographic directions 
(4 0), a total of 82% form a dimer before failure, while 
only 18% break from the monomer. In the cases studied 
in detail, the monomer is formed before the dimer, but, in 
a few cases, the contact breaks before forming the dimer. 



IV. TRANSPORT PROPERTIES 

We have finally computed the transport properties of 
the two possible types of configurations, i.e., a monomer 
and a dimer, before failure. The basics to calculate 
the zero-bias, zero-temperature conductance, G, in a 
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FIG. 6: (Left panel) Transmission for both spin species as a 
function of energy for the monomer configuration (shown in 
the inset) before failure. (Right panel) The same, but for a 
dimer configuration before failure. 

metallic nanocontact are contained in Landauer's formal- 
ism, where G is proportional to the quantum mechanical 
transmission probabihty of the electrons at the Fermi en- 
ergy, Ef 

G^^[T^{EF)+T^iEp)] (1) 

In this expression, the contributions from spin up (ma- 
jority) and spin down (minority) channels have been ex- 
plicitly separated, while the contribution from all the 
orbital channels has been condensed in T. For sim- 
plicity, we assume no spin mixing due to either spin- 
orbit scattering or noncoUinear magnetic structures at 
the bridge. The detailed electronic and magnetic struc- 
ture of the nanocontact is important, and, in order to 
achieve a quantitative level of agreement with experi- 
ments, one has to rely on first-principles or ab initio 
calculations. These calculations are performed with our 
code ALACANTiiii^ii^JI. The details of the calculation 
have been presented in previous publicationsi^ii^ii^iii. 
Essentially, one computes the self-consistent Kohn-Sham 



Hamiltonian for the narrowest part of the nanocontact, 
replacing the rest of atoms by a self-energy calculated 
based on a parametrized Bethe lattice. 

Using as input data two representative atomic configu- 
rations, as those shown in Fig. 6, the transmission spec- 
trum of these two structures has been calculated at the 
local spin density approximation (LSDA) level, and close 
attention has been paid to the choice of basis set for the 
central part of the nanocontact. As expected for Ni, ma- 
jority conduction is smooth as a function of energy due to 
the s-like nature of this channel, while minority conduc- 
tion is strongly fluctuating close to the Fermi energy due 
to the d-like character of this channeltf. Interestingly, 
the average value of the conductance around the Fermi 
level for both the examples lies somewhere in the vicinity 
of 1.6(2e^//i), which agrees fairly well with the value of 
the highest conductance peak in the histogram (see Fig. 
1). Remarkably, this value can only be obtained with an 
LSDA Kohn-Sham potential. The use of generalized gra- 
dient corrected functionals or hybrid functionals reduces 
strongly the conductance for minority electrons, leaving 
no possible explanation for the high-conductance peak. 

As far as the origin of the low-conductance peak, one 
could possibly attribute it to the presence of a domain 
wall at the narrowest section. Domain walls have a small 
but sizeable effect on the conductance of Ni nanocontacs, 
reducing it by an amount that agrees in magnitude with 
the conductance of the lowest peak. More experimental 
and theoretical work is, however, needed in this direction 
before this hypothesis can be confirmed. 
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